Solar Physics 
DOI: 10.1007/' 



Cross Helicity and Turbulent Magnetic Diffusivity in 
the Solar Convection Zone 



G. Riidiger^ • L.L. Kitchatinov^ • 
A. Brandenburg"''^ 

© Springer •••• 

Abstract In a density-stratified turbulent medium the cross helicity (u' • B') 
is considered as a result of the interaction of the velocity fluctuations and a 
large-scale magnetic field. By means of a quasilinear theory and by numerical 
simulations we find the cross helicity and the mean vertical magnetic field anti- 
correlated. In the high-conductivity limit the ratio of the helicity and the mean 
magnetic field equals the ratio of the magnetic eddy diffusivity and the (known) 
density scale height. The result can be used to predict that the cross helicity at 
the solar surface exceeds the value of 1 gauss km s~^. Its sign is anti-correlated 
with that of the radial mean magnetic field. Alternatively, we can use our result 
to determine the value of the turbulent magnetic diffusivity from observations 
of the cross helicity. 
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1. Introduction 

Dynamo theory for convective zones needs to know both the values of the a-effect 
and the eddy diffusivity. The a-effect is strongly related to the kinetic helicity 
which has opposite signs in the two hemispheres. Almost all of the theoretical cal- 
culations for rotating stratified turbulence lead to negative helicity {i.e. positive 
a-effect) for the northern hemisphere and positive helicity (i.e. negative a-effect) 
for the southern hemisphere. Despite all of the complications to measure the 
helicity on the solar surface, a new result has recently been presented by Komm, 
Hill, and Howe (2008). They do indeed find negative (positive) values for the 
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kinetic helicity in the northern (southern) hemisphere. This resuh is based on 
a ring-diagram analysis of GONG data. The values remained constant as long 
as the magnetic field did not exceed 10 gauss (G). Using the observation of the 
DIV-CURL correlation, C = ((m,^ + v^y)(v^x — u,y)), which is proportional to 
the kinetic helicity (Riidiger, Brandenburg, and Pipin, 1999; see Riidiger and 
HoUcrbach (2004) for more details), Duvall and Gizon (2000) found C to be 
negative (positive) in the northern (southern) hemisphere, as derived from the 
horizontal velocity components of mcsogranulation patterns. Egorov, Riidiger, 
and Ziegler (2004) simulated these observations with the NIRVANA code and 
reproduced them with Taylor numbers as small as 10'^. 

By use of the finding of Keinigs (1983) the current helicity and the a-efFect 
are anti-correlated, so one also can derive the sign of the a-effect by observation 
of the current helicity (J' • B'). Seehafer (1990) started to observe the current 
helicity at the solar surface showing that it is negative (positive) in the northern 
(southern) hemisphere. Again the a-cffect is found to be positive (negative) in 
the northern (southern) hemisphere. 

The numerical value of the helicity derived by Komm et al. (2008) is of the 
order of 10~^ cms~^. This turns out to be very small, because the resulting 
a-effect is below 1 cms^^. By comparison, Kapylii, Korpi, and Brandenburg 
(2009) find a ~ O.OSitrms near the surface from their convection simulations. 
With Wrms — 300 ms~^, this corresponds to 10 ms~^. The maximal a-value in 
their box center is of the order of O.Suims- This highlights a major discrepancy 
between theory and observations or, at least, a difficulty in determining a from 
observations. 

The empirical definition of the turbulent magnetic diffusion seems to be 
more straightforward. The decay of non-permanent magnetic structures such 
as sunspots or larger active regions lead to numerical values of the turbulent 
magnetic diffusivity. One finds Tyx — 10^^ cm^ from sunspot decay (Stix, 
1989) or rjT — 10^^ cm^ s~^ from the decay of active regions (Schrijver and 
Martin, 1990). These values are smaller than the value of 3 x 10^^ crn^s"^, 
which results from the widely used formula ryx ~ Mrms^/3 with correlation or 
mixing length £ and parameter values taken close to the surface. There is no 
possibility until now to observe the turbulent diffusivity on the solar surface for 
the quiet Sun where the magnetic quenching of this quantity by magnetic fields 
is negligible. We shall demonstrate in the present paper that there is a rather 
simple possibility to observe the magnetic diffusivity even in the presence of very 
weak magnetic fields (< 10 G), for which quenching should be negligible. 

We shall show that in the presence of a mean magnetic field along the direc- 
tion of density stratification, hydromagnctic turbulence will attain cross helicity, 
(u'-B'), whose value is proportional to the turbulent magnetic diffusivity. Indeed, 
our work is an extension of that by Kleeorin et al. (2003), who considered the 
effect of stratification of turbulent intensity. More recently, Kuzanyan, Pipin, 
and Zhang (2007) emphasized the importance of cross-helicity for estimating 
internal solar parameters important for the dynamo. We mention in this con- 
nection that cross helicity itself may constitute a potentially important dynamo 
effect dYoshizawa, 1990|[Yokoi, 1996[ ). In this paper we propose a measurement 
of cross-helicity in the Sun for estimating the turbulent magnetic diffusivity in 
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quiet regions. We argue that this is more accurate than measuring, for example, 
the mean electromotive force. 

2. Mean-Field Electrodynamics 

Let u' = u — (u) and B' = B — (B) be the fluctuations of velocity and mag- 
netic field about an average value denoted by angular brackets. The mean-field 
dynamo theory of cosmic magnetic fields is based on the relation 

(u' X B') = a(B) - /3(J) (1) 

between the turbulent electromotive force £ = {W x B') and the mean-field 
quantities (B) and (J), where J is the mean current density. For the purpose of 
this discussion we neglect here possible effects of mean flows on the correlators, 
i.e. we assume (u) = 0. Note the basic difference between the quantities a and (3 
in that /? is a scalar while a is a pseudoscalar. For rotating stars a pseudoscalar 
can be formed by use of the basic rotation rate, e.g., a oc (g • fZ) with gravity 
g as the only remaining preferred direction apart from H. Hence, the amplitude 
of the a-effect must mainly be influenced by the Coriolis number 

n* = 2Teorrf^, (2) 

where Tcorr is the correlation time of the dominating mode of turbulence. How- 
ever, J7* is very small at the solar surface so that the a-effect in Equation ([1]) 
cannot be observed directly. 

The parameter /? in Equation ([T|) exists even in nonrotating plasmas. It is thus 
not governed by the Coriolis number fl* and is therefore not small by comparison. 
It is, however, not possible to observe by direct means the mean current density 
(J) at the solar surface. The decay of sunspots should provide good estimates 
for /3 when the induction equation is solved by using Equation ([1} ; and the time- 
dependent solutions are compared with the observations (Krause and Riidiger, 
1975). Though successful, this procedure cannot serve as a proof of the existence 
of Equation ([ij. We must conclude, therefore, that the basic Equation ((T)) cannot 
be tested with observations taken from the solar surface. This is an unsatisfying 
situation given that Equation (II} is a fundamental relation of a whole branch of 
cosmic MHD and, of course, there is no better laboratory than the Sun to probe 
such basic relations. 

Fortunately, the situation is quite different for another correlation between 
fluctuations of flow and field, namely the cross hclicity (u' • B'), which itself is a 
pseudoscalar. It is straightforward to formulate the relation 

(u'-B') =ac(g-B)-/3e(n-J) (3) 

similar to Equation ([T|) . The ac-effect does not run with the Coriolis number n* . 
Similar to the a-effect in Equation ([T]) the ac in Equation ^ is of the dimension 
of a velocity but this velocity should be much larger than the corresponding a 
in Equation ((!]) . As the second term on the RHS of Equation ([3]) only exists in 
the presence of rotation, it will be negligibly small at the solar surface. 
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In summary, by simple reasons the observations of the cross correlation (u'-B') 
at the solar surface should give a realistic chance to confirm the existence of 
relations that are typical for mean-field electrodynamics. 



3. Nonconservation of Cross-Helicity in Turbulent Fluids 

The cross-helicity is conserved globally (as volume integral) in ideal incompress- 
ible fluids ( [Woltjer, 1958[ ) . In view of the conservation law, it may be anticipated 
that the balance of small-scale cross-helicity should be treated globally by defin- 
ing the small-scale cross-helicity sources, cross-helicity fluxes, and formulating 
the dynamical cross-helicity equation similar to the approach used to study the 
balance of magnetic helicity. Examples of such an approach to the cross-helicity 
problem can be found in the literature ( |Sur and Brandenburg, 2009p . In this 
section we show, however, that cross-helicity is not conserved in turbulent fluids 
such as the solar convection zone and its balance is controlled by local processes. 

The turbulence is known to dissipate efficiently the quantities, which are 
conserved in ideal fluids (with zero diffusivities). The well known example is the 
energy balance. Kinetic energy is conserved in ideal hydrodynamics. The rate 
of energy dissipation in Kolmogorov (1941) picture of turbulence is, however, 
constant independent of whatever small (but finite) is the viscosity. The same is 
true about almost all quantities conserved in ideal fiuids. Turbulent fragmenta- 
tion of scales cascades rapidly the quantities to the smallest scales where they 
dissipate. 

The only known exception is magnetic helicity that is conserved even in 
turbulent fiuids. The reason can be seen from the following. In the simplest 
case of isotropic homogeneous turbulence, the spectrum tensor of fluctuating 
magnetic fields can be written as 

B..(l,)~^fs.-'^]-i^e k (4) 
^^=^^>- 87rfc2 [^^^ fc2 ) 8^fc2 ^^J"'^"' ^4) 

where E^^^ and i?"^ are magnetic energy and helicity spectra, 

oo CO 

{B'^) = J E'^{k)dk, {B' ■ A') ^ J H'^{k)dk, B'^rotA'. (5) 



The spectrum tensor (j?]) is positive definite, BijCiC* > ( [Bochner, 1933[ ), where 
C is an arbitrary vector and the asterisk marks complex conjugation. For the 
tensor (j?]), this leads to the inequality ( [Moffatt, 1969 1 



\W{k)\k< E"'ik), (6) 

which is also known as the realizability condition. Imagine that at some (small) 
wavenumber ki helicity is finite and the spectral magnetic helicity to energy ratio 
at that wavenumber is k^^ = \H™{ki)\/E'^{ki). If helicity could follow magnetic 
energy in its cascade to large k, the ratio H^^/E™ would be constant across the 
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spectrum. Then, Equation ^ would require k < kg, whieli is an inequality 
that is impossible to satisfy for spectra with sufficiently broad inertial range. 
Therefore, magnetic helicity cannot be cascaded to diffusive scales to dissipate. 
The helicity is conserved, and the conservation law has important consequences 
for large-scale dynamos ( [Brandenburg and Subramanian, 2005[ ). 

The conservation of magnetic helicity is, however, an exception. For example, 
for kinetic helicity, which is also conserved in ideal hydrodynamics, instead of 
inequality ^ we have \H^\ < E^k with no restrictions for the kinetic helicity 
cascade to viscous scales. As a consequence, kinetic helicity is not conserved. The 
same is true for the cross helicity. This can also be seen by comparing the rates 
of hclicity dissipation. Using the fact that vorticity and current density scale 
inversely proportional to the square roots of viscosity and magnetic diffusivity, 
respectively, we see that the rate of magnetic helicity dissipation decreases with 
decreasing magnetic diffusivity proportional to its square root, while that of cross 
helicity is independent of viscosity and magnetic diffusivity and does therefore 
not vanish. The cross-hclicity balance is controlled by local processes. In spite 
of some striking similarities in the saturation of dynamos controlled by mag- 
netic and cross hclicity, the presence of significant cross-helicity dissipation as 
well as the forcing term in the momentum equation destroy the nice analogy 
( |Sur and Brandenburg, 2009[ ). In the following we proceed with deriving the 
cross-helicity from local relations. 

4. Quasilinear Theory of Cross Helicity 

In this section we derive the symmetric part {u[B'^Y = {WiB'j) + {u'jB[))/2 
of the cross correlation tensor {u'^B'^). The pseudotensor (u'^Bj) can be finite 
only in the presence of a mean magnetic field (B) and for inhomogcneous fluids. 
The required inhomogeneity can be due to stratification of density or turbulent 
intensity as well as the inhomogeneity of the mean field itself. 

The turbulent flow is assumed anelastic, so that div(pu') = 0. It is convenient 
to use the Fourier transformation of the momentum density m = pu', i.e. 



and similarly for the fluctuation of the magnetic field. The linearized equation 
for magnetic fluctuations in terms of the Fourier amplitudes reads 




(7) 



(-iw + r]k^)Bl{k,Lj) 




(8) 



where B is the Fourier transform of the mean magnetic field. 
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The spectral tensor of the momentum density that accounts for the stratifi- 
cation of the turbulence to first order terms reads 



X (% - hkj/e + {n,kj - Kjh) /{2k^)) , (9) 
where k = (z — z')/2, /« = z + z', g is the Fourier transform of the local spectrum, 

q{k,uj,r) ^ p^E{k,uj,r) ^ J q{k,u;, ny'" '' dn, (10) 

so that 



oo oo 

{u'^) = J J E{k,uj,r) dkduj. (11) 





Derivation of the cross correlation yields 

{u'.B'^r = iG.{B,) + G,{B,)) + (^^m + ^ry) 5., (U • (B)) + 
^r,T-^f/) mB,) + U,{B,)) 

+ ^7?) ((B,,.) + (B.,,)) , (12) 

where G = Vlogp and U = Vlog(u'^) are the gradients of density and turbulent 
intensity and 

oo oo ^ 



oo oo 

nk^uj^E , , 

' ;dkdoj. (14) 







(^2 + ^2;,4) 



Here rj = 1/ ^ou is the molecular magnetic diffusivity. Both quantities run with 
the magnetic Reynolds number for low conductivity (cr — > 0) and become finite 
for high conductivity (77 — >■ 0). From the cross correlation tensor the cross 
helicity 

(u' . B') = (G . (B)) + (^!^ + (U . (B)) (15) 

is obtained. 

Current observations only supply the correlation (it'^i?,',). From Equation (fT2|) 
we find 



= ,.G(S.) - ( ^ + i) - U{B^r) ) , (16) 
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where G = Gr and U = Ur are the only non-zero radial components of the strat- 
ification vectors. Further simplifications can be obtained by using the mixing- 
length approximation for the turbulence spectrum, 

E{k,Lu,r)=2{u'^)S{k-£-')S{uj), 7^ = f/r,,,, (17) 

(Kitchatinov, 1991), where £ is mixing length and Tcoit is the correlation time. 
It yields 

{<B'r) = VT (g{B,.) + ^U{B,.)^ . (18) 

The result can be explained as follows. A rising fiuid element > expands 
so that B'^. has the opposite sign as (B,.)- The fiuid particles which go down, 
uj, < 0, compress and B^ has the same sign as Br- The sign of the product u'^B'^, 
is opposite to (Br) in both cases - in accord with the first term on the right hand 
side (RHS) of Equation (|T5)) : note the negativity of G. An upward divergence of 
the mean field reduces the effect of density stratification. This is realized by the 
second term on the RHS of Equation ([T5| . The third term shows that also the 
non-uniformity of the turbulent intensity makes a contribution. However, the 
contribution of density stratification is dominant, because the density gradient 
in the upper convection zone is larger than the turbulent intensity gradient. 
This is already clear from Figure 1 of Krivodubskii and Schultz (1993), who 
plot, for a solar structure model, the relative contributions from G and U in 
the expression for the usual a effect, where both enter in equal amounts. Here, 
however, U enters with a 3/10 factor and is even more subdominant. Thus, we 
conclude that a finite cross correlation ([T5| indicates the presence of a large-scale 
radial field of the opposite sign. 

The leading term on the RHS of Equation (fT8)) is due to the density gradient. 
The resulting relation then reads 



{Br) H, 



p 



The magnetic eddy diffusivity can thus be determined if the LHS of Equation 
([TO)) is observed and the density scale height Hp is known from models of the 
solar atmosphere. 



5. Numerical Simulation 

It is straightforward to verify the validity of Equation (|19p using numerical 
simulations of isothermally stratified forced turbulence in a layer with constant 
gravity, g = (0, 0, —17) in Cartesian coordinates. In that case the scale height. 
Hp ~ c^/g, is constant. 

We perform simulations in a cubic domain of size , so the minimal wavenum- 
ber is fc = fci = 2it/L. We solve the governing equations of compressible mag- 
netohydrodynamics with an isothermal equation of state. The flow is driven by 
a random forcing function consisting of non-helical waves with wavenumbers 
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Figure 1. Dependence of the normalized cross helicity on Rm for various field strength 
B:,/Bcq < 0.1, Pm = 1, k{ / ki = 2.2, and Hpki = 2.5. The straight line denotes the fit 

(u'-B')/(u'-B')rcf = 0.05/?,n. 



whose modulus lies in a narrow band around an average wavenumber fcf (which 
corresponds to used in the previous section). We arrange the amplitude of 
the forcing function such that the RMS Mach number is around 0.1 or less, so 
the effects of compressibility are negligible. 

In all our runs we adopt stress-free pseudo-vacuum boundary conditions on 
the top and bottom boundaries, i.e. the horizontal magnetic field vanishes. The 
magnetic field is expressed in terms of the vector potential A as B = Bq-I- V x A, 
where Bg = (0, 0, Bq^) = const is the imposed vertical field which is fixed for each 
run. The simulations were performed with the Pencil Code0, which uses sixth- 
order explicit finite differences in space and third-order accurate time stepping 
method (Brandenburg and Dobler, 2002). A numerical resolution of up to 256^ 
meshpoints was used, depending on the value of the magnetic Reynolds number. 

We perform simulations for a number of different parameter combinations. 
The parameters that arc being varied include the strength of the imposed vertical 
field the forcing wavenumber fcf, the gravitational acceleration g, and hence 
Hp, and the values of the magnetic diffusivity. We express these quantities in 
non-dimensional form and define the magnetic Reynolds number as 



i?n 



r]k{ 



(20) 



^http : //pencil-code ■ googlecode ■ com | 
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Figure 2. Dependence 
strength for i?,m = 10, Pm 



of the normahzed cross hchcity on the normaUzed vertical field 
= 1, k{/ki = 2.2, and Hpk\ = 2.5. The dotted line corresponds to 



the graph of 0.3/[l + ((B)/Bi.ef )^]^''^ with S^ef = 0.85 Sc, 



The strength of the magnetic field is characterized by the mean equipartition 
field strength, 



-Soq = VJm){p} Urms, (21) 

which is of the order of 1000 G at the solar surface. We determine the cross 
helicity, (u' -B'), as a volume average. In order to relate this to Equation we 
also need to estimate the value of the turbulent magnetic difi^usivity. Earlier work 
showed that, to a good approximation, t^t can be estimated by Sur, Brandenburg, 
and Subramanian (2008), 

rjT ~ VTO ^ Urms/3fcf, (22) 

provided R,^ ^ 1, i.e. in the high-conductivity approximation. In a number of 
cases we have verified the validity of this approximation also for the stratified 
runs shown here. 

We present the value of (u' • B') in non-dimensional form by dividing by a 
reference value defined after Equation ([T^ as 



(u'-B'U^-^. (23) 

For small i?ni the normalized cross helicity depends on Rj^ (see Figure [T]) but 
it reaches unity for large i?ni- It is the expected behavior as Equation (j22p gives 
only a good approximation for Equation p3p for the case of high conductivity, 
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i.e. for cr — oo. In the opposite case of cr 0, expression p3p vanishes so that 
the small numbers of the lower-left corner of Figure [1] become understandable. 
For small values of i?in we have tjt oc R„^. For the largest values of i?ni the error 
bars for the numerical results are larger. This is mainly because those simulations 
require larger numerical resolution and long run times become prohibitive. 

Figure [2] shows the dependence of the normalized cross helicity, defined by 
the ratio (u' • B')/(u' • B')i.of, on the normalized field strength Bo^/Bcq. Note 
that the cross helicity is quenched by nearly a factor of 10 for Bqz ~ Bcq. 



6. Conclusions 

We have shown that nonrotating turbulence at the top of the solar convection 
zone under the influence of a vertical magnetic field yields a finite cross helicity. 
The only requirement is the existence of density stratification which enters the 
induction equation via the anelastic relation div(/ou') = 0. The Boussinesq ap- 
proximation cannot be used. The effect exists mainly in the high-conductivity 
limit, i.e. for sufficiently large magnetic Reynolds numbers (see Figure [T]). The 
radial magnetic field, on the other hand, must be weak enough to remain passive 
so that it does not dominate the flow. Figure [5] shows that the maximum field is 
given by -Bgq which is much higher than the mean vertical field, which is of the 
order of a few gauss on the Sun. 

To estimate the value of the cross helicity at the solar surface we shall assume 
a density scale height of 100 km. Then one finds from Equation (fTO|) that 

«S;>--^|jGkms-^ (24) 

The average is here to be applied over many turbulent cells which, in the Sun, 
might correspond to 30-100 Mm. The magnetic diffusivity in Equation has 
been used in the form ryx = 10^^7712 cm^ s~^ and the density scale height as 
Hp = lOOi/y km. We thus predict the existence of a cross helicity of more than 
1 Gkms^^. We also emphasize that the cross helicity is anti-correlated to the 
mean radial magnetic field, i.e. 

{u[.B[.) {Br) < 0. (25) 

For a dipolar background field the sign of the cross helicity will be opposite in 
the two hemispheres. 

Relation can also be used to measure the magnetic diffusivity if the cross 
helicity is known by observations. In order to find the cross helicity one only has 
to correlate the observed flow fluctuations with observed magnetic fluctuations. 
Together with the calculated mean value of the radial magnetic field. Equation 
(HH) provides the unknown quantity rjT- We hope that such an analysis of the 
observations using, for example, data from the Hinode satellite will soon provide 
supporting evidence for an anti-correlation between (m^S^) and (Br), and that 
a meaningful value of rj^ can be obtained in that way. 
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